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Abstract 

We consider the regularity of the extremal solution of the nonlinear 
eigenvalue problem 

( >X \ u = on on, 

where O is a smooth bounded domain in and c{x) is a smooth bounded 
vector field on Q. We show that, just like in the advection-free model 
(c = 0), all semi-stable solutions are smooth if (and only if) the dimen- 
sion N < 7. The novelty here comes from the lack of a suitable variational 
characterization for the semi-stability assumption. We overcome this dif- 
ficulty by using a general version of Hardy's inequality. In a forthcoming 
paper [5], we indicate how this method applies to many other nonlinear 
eigenvalue problems involving advection (including the Gelfand problem), 
showing that they all essentially have the same critical dimension as their 
advection-free counterparts. 

1 Introduction 

The following equation has often been used to model a simple Micro- Electro - 
Mechanical System (MEMS) device: 

, , f -Am = n A V2 in fi, 

(P)x < { u) 

1 M \ u = ondil, 

where 17 is a smooth bounded domain in M. N , A > is proportional to the 
applied voltage and < u(x) < 1 denotes the deflection of the membrane. 
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This model has been extensively studied, see [TU], [H] in regards to the model 
and [7], [6], [8] for mathematical aspects of (P)x- It is well known (see above 
references) that there exists some positive finite critical parameter A* such that 
for all < A < A*, the equation (P)x has a smooth minimal stable (see below) 
solution u\, while for A > A* there are no weak solutions of {P)\ (see [7] for 
a precise definition of weak solution) . Standard elliptic regularity theory yields 
that a solution u of (P)x is smooth if and only if sup n u < 1. One can also show 
that A i— > u\(x) is increasing and hence one can define the extremal solution 

u*(x) := lim u\(x), 

which can be shown to be a weak solution of (P)\* . 

Recall that a smooth solution u of (P)\ is said to be minimal if any other 
solution v of (P)a satisfies u < v a.e. in Q. Such solutions are then semi-stable 
meaning that the principal eigenvalue of the linearized operator 

L UiX := -A - 



(1-u) 3 

in Hq(Q,) is nonnegative. This property can be expressed variationally by the 
inequality 

2A / 7-^3 < / |vvf, W> e ir^fi). (i) 
Ai (i - wr Jo 

which can be viewed as the nonnegativeness of the second variation of the energy 
functional associated with (P)\ at u. 

Now a question of interest is whether u* is a smooth solution of (P)\* ■ It is 
shown in [7j that this is indeed the case provided N < 7. This result is optimal in 
the sense that u* is singular in dimension N > 8 with Q taken to be the unit ball. 

Our main interest here will be in the regularity of the extremal solution 
associated with 



-Au + c{x) ■ Vu = (i_ u \a m 

u = on dfl, 



where c G C°°(il,R N ) and where again is a smooth bounded domain in M. . 
Modifying the proofs used in analyzing (P)\ one can again show the existence 
of a positive finite critical parameter A* such that for < A < A* there exists 
a smooth minimal solution u\ of (S)\, while there are no smooth solutions of 
(S)\ for A > A*. Moreover, the minimal solutions are also semi-stable in the 
sense that the principal eigenvalue of the corresponding linearized operator 

2A 

i«,A,c := —A + c(x) ■ V — 



(1 - ux) 3 



in i?o(fi) is non-negative. See [3] where these results are proved for general C 1 
convex nonlincarities which are super linear at oo. Our main result concerns the 
regularity of the extremal solution of {S)x- 
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Theorem 1. If 1 < N < 7, then the extremal solution u* of (S) a* is smooth. 

Remark 1. A crucial (in fact the main) ingredient in proving the regularity of 
u* in (-P)a, is the energy inequality |T]) which is used in conjunction with the 
equation (P)\, to obtain uniform (in A) L p -estimates on (1 — u\)~ 2 whenever 
u\ is the minimal solution (See [7]). However, the semi-stability of u\ in the 
case of (S)\, does not translate into an energy inequality which allows the use 
of arbitrary test functions. Overcoming this will be the major hurdle in proving 
Theorem Q] 

We point out, however, that if c(x) = V7 for some smooth function 7 on f2, 
then the semi-stability condition on the minimal solution u\ of (S)\ translates 
into 

2A / 7T— ¥vi ^ I e ~ 7 l v ^l 2 > W> e H^Sl). (2) 
Jn K 1 — u \) Jn 

Then, with slight modifications, one can use the standard approach for (P)\ to 
obtain the analogous result for (S)\ stated in Theorem [TJ 

The novel case is therefore when c is a divergence free vector field. Actually, 
we shall use the following version of the Hodge decomposition, in order to deal 
with general vector fields c. 

Lemma 1. Any vector field c S C 00 ^,!^) can be decomposed as c(x) — 
— V7 + a(x) where 7 is a smooth scalar function and a(x) is a smooth bounded 
vector field such that div(e 7 a) = 0. 

Proof. By the Krcin-Rutman theory, the linear eigenvalue problem 

J Aa + div(ac) = /la Q, , . 

\ (Va + ac)-n = dfl, ^ 

where n is the unit outer normal on d£l, has a positive solution a in Q when /j, is 
the principal eigenvalue. Integrating the equation over f2, one sees that /i = 0. 
The positivity of a on the boundary follows from the boundary condition and 
the maximum principle. In other words, we have that Aa + div(ac) = on fi, 
and a > on f2. 

Now define 7 :— log(a) and a := c + V7. An easy computation shows that 
div{e~<a) =0. □ 

Throughout the rest of this note c, a, 7 will be defined as above. 



2 A general Hardy inequality and non-selfadjoint 
eigenvalue problems 

Consider the linear eigenvalue problem 

r -A(j> + c ■ - 9 <t> = Kef, n, 

\ cj) = an, [ > 
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where c is a smooth bounded vector field on Q, p € C°°(£l) and K is a scalar. 
We assume that (4>,K) is the principal eigenpair for ((4]) and that <fi > in fi, 
and if > 0. Note that elliptic regularity theory shows that <f> is then smooth. 

We shall now use a general Hardy inequality to make up for the lack of a vari- 
ational characterization for the pair (<fi, K) . The following result is taken from 
[4 , which we duplicate here for the convenience of the reader. For a complete 
discussion on general Hardy inequalities including best constants, attainability 
and improvements of, see [3] . We should point out that this approach to Hardy 
inequalities is not new, but it is generally restricted to specific functions E which 
yield known versions of Hardy inequalities; see [T] and reference within. 

Lemma 2. Let A(x) denote a uniformly positive definite N x N matrix with 
smooth coefficients defined on fl. Suppose E is a smooth positive function on Q 
and fix a constant (3 with 1 < f3 < 2. Then, for all ip G ^o(^) we have 

J Si J O 

where J Q |V^|i = J n A{x)Vip ■ Vip. 

Proof. For simplicity we prove the case where A(x) is given by the identity 
matrix. For the general case, we refer to [4]. Let Eo denote a smooth positive 
function defined in fl and let ip € C* c °°(ft). Set v := -f=. Then 



|V£r 12 
AE? 



|V?/f =E \Vv\ 2 + ^-^-^ 2 +vVv VE . (6) 



Integrating the last term by parts gives 



/ 

Jn 



vVWE = 1 -[ 
n z Jn 



and so integrating (O gives 

|V^>i/^ + i/^V, (T) 
n 4 Jn A) 1 Jn 

where we dropped a nonnegative term. So we have the desired result for (3 = 1. 
When [3 ^ 1 one puts Eo :— E 13 into (O and collects like terms to obtain the 
desired result. □ 

We now use the above lemma to obtain an energy inequality valid for the 
principal eigenpair of Q. 

Theorem 2. Suppose that the principal eigenpair (</>, K) of ^ are such that 
4> > and K > 0. Then, fori < f3 < 2 we have for all ip £ H^[n), 

' . w > s?f£ / I [ *^V. is) 



<> 4 J n 4> 2 2 J n 2 Jn 
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Proof. Note that ^ can be rewritten as 

-div(e 7 V0) + e 7 a • V<j> = e 7 (p(x) +K)(f> in 0, 

where as mentioned above we are using the decomposition c = — V7 + a. We 
now set E := <j) and A(x) — eT (where / is the identity matrix) and use ([5]) 
along with the above equation to obtain the desired result. Note that we have 
dropped the nonnegative term involving K. 

□ 



3 Proof of theorem [T] 

For < A < A*, we denote by u\ the smooth minimal semi-stable solution of 
(S)\. Let (4>, K) denote the principal cigenpair associated with the linearization 
of {S) x at u\. Then < 4> in 0, < K and (4>, K) satisfy 

Again, elliptic regularity theory shows that <f> is smooth. Consider c = — V7 + a 
to be the decomposition of c described in Lemma 1. We now obtain the main 
estimate. 

Theorem 3. For < A < A* ; 1 < (3 < 2 and < t < (3 + y/ f3 2 + (3, we have 
the following estimate: 



X+lJ Ja(l-n)*+ 3 - in (1 - «a)'« H2-flJaQ- «>.)*' 

Proof. Fix 0</3<2, letO<t and u denote the minimal solution associated 
with (S)x- We shall use Theorem H with p(x) = n^-p- Put ip := - 1 

into (jHJ to obtain 

^ / ? e7 ' v t 2 +9 > /3A [ ( _l_ _ r 2 

i/ e1 / (2-/3)|V^| 2 a-Vfy 2 
" n V 2 2 



2 

Now note that (S)\ can be rewritten as 

Ae 7 

— div(e 7 Vu) + e 7 a • V« = — in fi, 

(1 — v,y 

and test this on cf> := ( - 1 _ tt 1 ^ t +i — 1 to obtain 

{2t + 1) L a - «)"' + 2 + H = A L o^F ( Jy=W- 
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where 

H := I e 7 a • Vw I - - 1 



o v(l-«) 2 
One easily sees that H = after considering the fact H can be rewritten in the 
form J n (e 7 a) • VG(u) for an appropriately chosen function G with G(0) = 0. 
Combining the above two inequalities and dropping some positive terms gives 



2t + lJJ n (l- u ) 2t + 3 " " J n (i- u )a+t 

P f IK, J 1 



e<k(x) - - 1 

2j n { ' {(1-uY 
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where 



Simple calculus shows that 

sup A (a;) < 



n 



2(2-/3)' 



which, after substituting into the above inequality, completes the proof of the 
main estimate. □ 

Note now that the restriction t < (3 + \J ' f3 2 + (3 is needed to ensure that the 
coefficient j3— is positive. It follows then that ^_ ux yi is uniformly bounded 

(in A) in L P (Q) for all p < po := | + y/6 « 5.94... and after passing to limits we 
have the same result for the extremal solution u*. 

To conclude the proof of Theorem [TJ it suffices to note the following result. 

Lemma 3. Suppose 3 < N < 7 and the extremal solution u* satisfies jr^ryi G 
L^(O), then u* is smooth. 

Proof. First note that by elliptic regularity one has u* G W 2 -^(fl) and after 
applying the Sobolev embedding theorem one has u* G C°^{Q). Now suppose 
— 1 so that there is some xo G fi such that u(xq) = 1. Then 

1 C 

> 



l-u{x) |a;-a:o|f 
and hence 

oo > / —^r > C 



a {{i-u*yy-t- ~ J a \x 



N 



which is a contradiction. It follows that ^_ u ,y <= L°°(Q), and u* is therefore 
smooth. □ 
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Using this lemma and the above L p -bound on , one sees that u* is 

smooth for 3 < N < 7. To show the result in dimensions N = 1,2, one needs a 
slight variation of the above argument. We omit the details, and the interested 
reader can consult [7] for the proof when c{x) = 0. 

Remark 2. As mentioned in the abstract, this method applies to most non- 
selfadjoint eigenvalue problems of the form 

J -Am + c(s) • Vw = Xf(u) in O, 
[ )x \ u = on <90, 

where f(u) is an appropriate convex nonlinearity such as f(u) = e u , and f(u) = 
(l+u) p . It shows in particular that the presence of an advection does not change 
the critical dimension of the problem, hence addressing an issue raised recently 
by Berestycki et al [3]. One can also extend the general regularity results of 
Nedev [9] (for general convex / in dimensions 2 and 3) and those of Cabre and 
Capella [2] (for general radially symmetric / on a ball, and up to dimension 9) . 
All these questions are the subject of a forthcoming paper [5]. 
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